submitted to acta physica slovaca 



1-6 



COHERENT PHOTOPRODUCTION OF 7/^ AND HIGH-MASS e+e PAIRS 
IN ULTRA-PERIPHERAL AU+AU COLLISIONS AT ^/TZ = 200 GeV 



David d'Enterria 
for the PHENIX Collaboration 

Nevis Labs, Columbia University, Irvington, NY 10533, and New York, NY 10027, USA 
LLR, Ecole Poly technique, 91128 Palaiseau, France, EU 

Submitted Dec 30, 2005 

Higii energy ultra-peripheral collisions (UPC) of heavy-ions generate strong electromagnetic 
fields which open the possibility to study YY and Y-nucleus processes in a kinematic regime 
so far unexplored. We report on preliminary PHENIX results of J/\\t and high-mass e+e~ 
photoproduction at mid-rapidity in coherent electromagnetic Au-l-Au interactions at ,ys^ = 
200 GeV tagged with forward neutron emission from the Aii* dissociation. 

PACS; 13.40.-f, 13.60.-r, 24.85.-l-p, 25.20.-x, 25.20.Lj, 25.75.-q 

1 INTRODUCTION 

1.1 YY, jp physics at e+e and ep colliders 

Our knowledge about elementary particles and their fundamental interactions is mainly obtained 
through the study of high energy particle collisions. A part from the more "conventional" e+e^, 
ep (DIS) and pp, pp collisions, high-energy YY 1 1 1 and yp [2 \ processes provide an interesting 
and complementary approach to study strong and electro-weak interactions. At high energies, a 
(quasi) real photon can interact as a point-like particle (e.g. Compton-like scattering) or through 
quantum fluctuations into a charged-fermion pair {qq, l^l^), a charged-boson pair (e.g. W+W^) 
or directly into a vector meson (carrying the same quantum numbers j'^^ = 1^ ^ as the photon). 
At energies below the VK+W^ threshold, the photon wave function can be written as \2i 

Iy)=co|Yo) + E ^y\^)+ L c,\qq) + c,\l+l-). (1) 

V=p,(0.(|),i/y.T q=u.cl,s.c,b l=e,/.i,z 

Though in practice the first term dominates (i.e. cq « 1), hadronic/partonic fluctuations (under 
the non-perturbative "Vector-Dominance-Model" form Y — > V or, at larger qq virtualities, via "re- 
solved" quark pairs qq which can, at their turn, radiate also gluons) interact strongly and give 
the largest contribution to the total YY- IP cross-sections at high energies. Fig.^shows the total 
hadronic cross-sections measured in e+e^, YY ™d jp collisions as a function of the center-of- 
mass energy y/s. While, the (s-channel) e+e^ annihilation cross-section decreases with increas- 
ing energy (Oee^hadnms °^ OCgm/'^)' photon-induced collisions behave like pp, pp collisions. Their 
cross-sections rise monotonically with ^/i in agreement with an increasingly large Pomeron (or 
two-gluon) exchange contribution as described by Regge phenomenology |3lEl- Although in 
the present-day there are no beams of high energy real photons available, beams of high energy 
quasi-real or virtual photons are commonly generated at e+e^ and ep colliders. Indeed, the 
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Fig. 1. Total hadronic cross-section me^e , yyand yp collisions as a function of the center-of-mass energy 
•Js. The data is from the PDG O, and the fits to yy, yp data are the latest Regge-theory parametrizations 
of the COMPETE collaboration |4|. 

electromagnetic field of any relativistic charged particle can be described as a flux of "equiva- 
lent" photons with energy distribution fy/eiz) — dNy/d(0 given by the Weizsacker- Williams |6| 
or Equivalent Photon Approximation (EPA) formula, which for an beam reads: 

^-^-[l + (l-^)']ln§^ ,with z = co/£.. (2) 

Here z is the fraction of the energy carried by the photon and = —{q^) the y virtuality. 
Qmin ~ "^e2^/(l — z) ~ GcV^/c^, whcfcas QI,^,^ depends on the properties of the produced X, 
e.g. Qf„^,^- ~ nip for hadron production, and « forX ~ l^l^ , where — Wyy = 4c0iC02 
(or ~ Wyp = AEpEez) is the squared yy (or yp) cm. energy. Fig. |2] shows the typical two- 
photon (e+e^ e+e^yy— + e^e^X) and photon-proton (ep ejp epX) processes commonly 
studied at LEP and HERA colliders. Two basic properties characterize yy and jp collisions: 

• The variable-energy photon "beam" always carries less energy than the parent lepton beam 
(the spectrum (|2j has a soft bremsstrahlung shape, dNy/dd) °^ ln{l / 03) / (a, peaked at low 
CO) and thus the available center-of-mass energies are typically 0.1 -y/s < Wyy_yp < 0.5\/s. 

• The (almost) real photons beams have, by definition, « (they are often kinematically 
selected by requiring the scattered parent to be close to the incident direction) and thus 
the produced X have always small transverse momentum (at variance with DIS scattering 
where the y* has 2^ > 2 GeV^/c^ and the produced particles have large pr). 

The main interest of y-induced processes at e^e^ ,ep colliders so far is that they allow for preci- 
sion studies of QCD dynamics (tests of Regge theory, low-energy quark model spectroscopy, 
BFKL evolution, hard-scattering factorization, diffractive interactions, ...) in a background- 
free environment and with a comparatively simpler initial state than in hadron-hadron colli- 
sions. QED processes (e.g. /+/^ production), on the other hand, have been less addressed 
than QCD-related studies. Typical physics measurements encompass |7|: (i) total hadronic 
cross-sections, (ii) yy widths (Tyy) of C-even (scalar 0^+'++ and tensor 2++^- - ) mesons, (iii) 




Fig. 2. Typical diagrams for photon-photon (left) and photon-proton (right) reactions at e e and ep 
colliders respectively. 

double vector meson production (p,a),(|),7/\(/,T), (iv) structure function of the photon, (v) (hard) 
photo-production of quarkonia (y/\(/,T), open heavy quarks and jets; and (vi) diffractive structure 
functions. The corresponding two-photon {<3yy^x) or photo-production {Oyp^x) cross -sections 
are theoretically computed using different methods depending on the cm. energy W and the 
nature of the produced particle X. Regge-based models (as e.g. in the DPM |8|), (General- 
ized) Vector-Meson-Dominance (VDM) approaches (51, pQCD factorization 1101 (with different 
available parametrizations of the parton distribution function of the photon fTT|), color-dipole 
(or saturation) approaches [12 131 . or a combination (of some) of them (as e.g. in PHOJET L14J 
or PYTHIA jj] event generators) are often used. 

At future linear collider facilities using beams of back-scattered laser photons off beam electrons, 
new prospects appear both within the Standard Model (Higgs production, triple gauge couplings, 
precision top physics) and beyond-SM (SUSY particle production, extra-dimensions, ...) F15I. 

1.2 yy, yA physics in ultra-peripheral Ah- A reactions 

As aforementioned, any relativistic charged particle generates a flux of quasi-real photons which 
can be used for photoproduction studies. Though in the previous Section we have considered the 
case of beams, the same holds obviously true for accelerated protons 1 1 6 1 or nuclei 1 17.18^19J . 
The corresponding Equivalent Photon Approximation formula for an extended object such as 
a proton or a nucleus with charge Z and mass ;«a is 12 

/y/AW = ^ + r \F{Q')?dQ\ (3) 

where is the momentum transfer from the projectile, Q^^j^ = m^z^/{l — z), and F{Q^) a nu- 
clear form factor describing its spatial distribution. From Eq. Q one can see that the equivalent 
photon flux in relativistic heavy-ion collisions (with charge Z~ 80 for conventional Au or Pb 
beams) has an amplification factor of ^ 6.5-10^ compared to or proton beams. The corre- 
sponding yy cross-sections are a factor Z^ ~ 4-10^ larger ! This fact drives the main interest on 
coherent electromagnetic (or ultra-peripheral) interactions of heavy-ions ll7lll8l[T9l . For colli- 
sions between heavy-ions, rather than Eq. Q it is more useful the expression for the equivalent 
photon spectrum above a given minimum impact parameter 1 20 1 which is usually chosen as twice 
the nuclear radius, b > fe„„„ k,2Ra (i.e. /7„„„ w 15 fm for Au, Pb beams), so that the contribution 
from non-electromagnetic ion-ion interactions can be removed: 

fy,A{z,b > b„u„) = ^^£!!1^1 [2xKo{x)Ki{x)~x^ {kI{x) - K^{x))] , (4) 
7Z z 
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where Ko i are the modified Bessel functions of the 2nd kind and x = zniAbmm- The vaHdity of 
the apphcation of the EPA formula for heavy-ions is Hmited to the case where all the protons of 
the nucleus interact electromagnetically in a coherent way. In that case, the wavelength of the 
resulting photon is larger than the size of the nucleus, given by its radius Ra- This "coherence" 
condition limits the maximum virtuality of the produced photon to very low values L17J : 

= (co^/T^ + ?i) ^ 1 /^A (where y is the beam Lorentz factor), (5) 

and thus for most purposes (but for lepton pair production) those photons can be considered as 
(quasi) real with maximum energy and perpendicular momentum: 

Y ^ 1 

« < (a„ax , and < ~ 30 MeV. (6) 

At RHIC (LHC) energies with y = 100 (y « 2800 for Pb), the maximum photon energy in the 
lab system is a)„K„ « 3 GeV (80 GeV). Thus, the corresponding maximum energies in the cm. 
system for yy and jA processes are Wyy"'^ « 6 (160) GeV and Wy]^^ « 34 (940) GeV respectively. 

1.3 Production cross-sections in UPC collisions 

The production cross-section of a system X in UPC Ah-A collisions is computed as the convo- 
lution of the corresponding equivalent photon spectrum @ with the elementary photonuclear or 
two-photon cross-sections: 

a{A+A^y+A+A^A+A+X) = J^^ fy/A{z)ciyA^x {WyA)dz (7) 

a{A+A^yy+A+A^ A+A+X) ^ fy/A{zi)fy/A{z2)Oyy-^x{Wyy)dzidz2i^) 

where, as aforementioned, the elementary OyA^x and Oyy^x are theoretically computed using 
different methods depending on the cm. energy W and the particle X. The two cases of interest 
in this paper, shown in Figure|3 are the coherent photoproduction of: 

I. //v|/, the heaviest vector meson effectively accessible in yA collisions at RHIC, via: 

A + A y+ A) A* + A W 

II. High mass di-electron continuum in yy collisions: A +A y+y) A* +A^*^ +e+ + e^. 

For the first process (I), OyA^y/yA has been determined using perturbative QCD for the jp 
J/'^p process ^21] plus initial-state shadowing and final-state 7/\|/-nucleus interaction [22], or 
within the color-dipole formalism |13j|. The cross-section for heavy vector meson (J/\\t,T) pho- 
toproduction is found to depend (i) quadratically on the gluon density Ga{x,Q^) | 23 | : 

- -^^16713 [xGA{x,Q^)f , with /4 , and x -M^/W.^^, (9) 

as well as on (ii) the probability of rescattering or absorption of the QQ pair as it traverses the 
nucleus. The study of quarkonia production in yA collisions at RHIC or LHC energies is, thus, 
considered an excellent probe of (i) the gluon distribution function Ga{xtQ^) in nuclei, and (ii) 
vector-meson dynamics in nuclear matter 

For the second purely electromagnetic process (II), the total cross-section is huge 

(about 30 kb for Auh-Au at RHIC !) but strongly peaked at forward-backward rapidities and 
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Fig. 3. Lowest order Feynman diagrams for y/\|f (left) and dielectron pair (right) production in YY and yA 
processes accompanied by Au Coulomb excitation in ultra-peripheral Au+Au collisions. 

at low m,„,. < 10 MeV/c^ [171. Several calculations exist for multiple pair production within 
QED 1 17 1 or solving the semi-classical Dirac equation |24|. The main interest of such coherent yy 
collisions is that one can test QED in a very strong field regime {Z(Xem ~ 0.6) where perturbative 
calculations are expected to break |25|. 



At RHIC energies, the production cross section for vector mesons (p, (|), T) in UPC AuH-Au 
is large and accounts for as much as ~10% of the total CauAu ~ 7 b nuclear cross section I26I27I . 
Measurements of coherent photonuclear production of the p meson |'28l as well as yy production 
of low-mass pairs L29J have been performed by the STAR collaboration. The PHENIX anal- 
ysis presented here aims at similar measurements but at higher invariant masses. 

Detailed knowledge of the experimental signatures of coherent UPC events is a basic pre-requisite 
to setup an efficient UPC trigger and to define the reconstruction and analysis cuts used in the 
present work. The typical characteristics of UPC AuH-Au events are: 

1. Low central multiplicities: typical values are (well) below ^15 particles. 

2. Low total transverse momentum ("coherence condition"): pj < 2hc/R w 50 MeV/c or 
Pt ^ minvh~ 30 MeV/c. 

3. Large probability of multiple electromagnetic interactions (y-exchanges shown in Fig.|3} 
leading to (single or double) nucleus Giant-Dipole-Resonance (GDR) excitation followed 
by Xn neutron(s) decay. Typical probabilities are P\„ ^30-50% (7/\|/) or Pi„ ^^20% (p) 
which factorize out when determining the UPC cross-sections I25II30I . 

4. Zero net charge: even number of charged tracks of opposite signs. 

5. Narrow dN/dy peaked at mid-rapidity (narrower for larger m,,,,,). 
Property 3. is the most useful for UPC tagging and trigger purposes. 



The data presented here were collected with the PHENIX detector at BNL RHIC during the 2004 
high luminosity AuH-Au run at ^/s^ = 200 GeV (Run-4). The gold beams had 45 or 56 bunches 
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with 10 ions/bunch and 106 ns crossing-time. Typical luminosities at the beginning of the store 
were ~2 lO^'' cm^^s^\ twice larger than the design luminosity (thanks mainly to an improved 
vacuum system) 1311 . The PHENIX detector 1321 is specifically designed to measure hard probes 
by combining good mass and particle identification (PID) resolution in a broad momentum range, 
high rate capability, and small granularity. The central-arm detectors used in this analysis (DC, 
PC, RICH and EMCal) are those needed for the measurement of J/\\t (e+e^ decay mode) and 
high-mass (m,„y > 2 GeV/c^) dielectron pairs at y = 0. The ultra-peripheral AuH-Au events were 
tagged by neutron detection at small forward angles in the Zero-Degree-Calorimeters (ZDC). 

The momentum and trajectory of the tracks were reconstructed with the central tracking (CNT) 
system (covering Ar| = 0.7 and A(|) = n) based on a multi-layer drift chamber (DC) followed by 
multi-wire proportional chambers (PC) with pixel-pad readout, both placed in an axial magnetic 
field parallel to the beam {B„,ax = 115 T m). Electron-positron identification was done with the 
Ring-Imaging- Cerenkov (RICH, with CO2 gas radiator) and electromagnetic calorimeter EM- 
Cal |331, with a total total solid angle at mid-rapidity of Ar| w 0.7 and A(|) = n. The PHENIX 
EMCal consists of six sectors of lead-scintillator sandwich calorimeter (PbSc, 15552 individ- 
ual towers with 5.54 cm x 5.54 cm x 37.5 cm, I8X0) and two sectors of lead-glass Cerenkov 
calorimeter (PbGl, 9216 modules with 4 cm x 4 cm x 40 cm, 14.4Xo), at a radial distance of 
^5 m from the beam line. The ZDC |34''35l hadronic calorimeters placed 18 m up- and down- 
stream of the interaction point, cover the very-forward region, |9| < 2 mrad, and measure the 
energy of the neutrons coming from the Au* Coulomb dissociation, with ^^20% resolution. 

The total equivalent sampled luminosity for this study was: i;,,, = Nbbc-lli /<^AuAu x £bbc — 
120 ± 10 /jb^\ where Nbbc-lli = 886 • 10^ is the total number of events collected (after vertex 
cut and QA) by the Beam-Beam-Counter (BBC) minimum bias AuH-Au trigger, and we use for 
this preliminary analysis the "nominal" values for the AuH-Au total cross-section {<3auAu = 6.85 
b) and BBC-LLl ti'igger efficiency {Ebbc = 92 ± 3%) l36l . 

2.2 IVigger and data sample 

The events used in this analysis were collected by a level-1 UltraPeripheral (LLl-UPC) trigger 
set up for the first time in Run-4, with the following characteristics: 

1. Veto on coincident signals in both BBC (covering 3.0 < |r|| < 3.9 and full azimuth) is 
imposed in order to reduce peripheral nuclear and beam-gas collisions. 

2. A large energy {E > 0.8 GeV) cluster in EMCal is required to select the from the J/\\t 
decay (E « 0.5mj/y « 1.6 GeV) and from the high-mass di-electron continuum. 

3. At least 30 GeV energy deposited in one or both of the ZDCs is required to select AuH-Au 
events with forward neutron emission (Xn) from the (single or double) Au* decay. 

The trigger used to detect high-energy electrons in the central arm was the standard EMCal- 
RICH-Trigger (ERT) with 2x2 tile threshold at 0.8 GeV. A software algorithm performed a 
crude reconstruction of EMCal clusters by summing the pedestal-subtracted and gain-calibrated 
energies in overlapping "tiles" of 2 x 2 calorimeter towers with an estimated (preliminary) effi- 
ciency for high-energy e"^ coming from the //\|/ decay of ejf^'^J^, — 0.9 ± 0. 1 . 

The total number of events collected by the UPC trigger was 8.5 M (i.e. 0.4% of the min.bias 
AuH-Au trigger) and the data set comprised ~ 1000 raw-data files (~ 1 . TB). Standard QA further 
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reduced this number down to 6.7 M events. Most of these UPC-triggered events were not, how- 
ever, signals of interest in this analysis {coherent J /\\t and high ). Other events (with 
their discriminating characteristics indicated) were likely to fire also the UPC trigger I17I18I19I : 

(1) Cosmic rays: no ZDC, no good vertex. 

(2) Beam-gas collisions: no good vertex, large multiplicity, asymmetric dN/dy. 

(3) Peripheral nuclear collisions: comparatively large particle multiplicities, large pair pr- 

(4) (Hard) diffractive collisions: Pomeron-Pomeron events with rapidity gap(s) f37"381 are 

usually accompanied by forward proton emission, pt{IPP) > PriyY), and like-sign pairs. 

(5) Coherent yy qq: produce mainly hadrons (removable by the identification cuts). 

(6) Incoherent (or "quasi-elastic") photon-nucleon (yA^) collisions: prijIP) > Priyj), wider 

and asymmetric dN/dy, Xn > 2n (recoiling nucleon induces nuclear break-up with larger 
probability than in fully coherent interactions) 1221 . 

Processes (1) and (2) can be considered "non-physical" sources of background and are rejected 
by simple global event analysis cuts. "Physical" processes (3), (4), (5) can be removed by offline 
cuts more easily than the yA^ contribution (6) which has experimental signatures very similar to 
true yA reactions (see 1221 and Section l33t . 

2.3 Data analysis 

Charged particle tracking in PHENIX central arms is based on a combinatorial Hough transform 
in the track bend plane, with the polar angle determined by PCI and the collision vertex along the 
beam direction 1 39 1. The original z- vertex of the track was measured with 1-cm resolution using 
a specific method based on PC hits and EMCal clusters since the standard procedure based on 
BBC and ZDC is not applicable for UPC events which, by definition, do not have BBC coinci- 
dences and often do not have ZDC coincidences. Track momenta are measured with a resolution 
bp/p w QJ%+l.Q%p w 1.7 - 2.7% for the relevant range (p ^ I - 3 GeV/c) in this analysis. 

The following global cuts were applied to enhance the sample of genuine y-induced events: 

1. A standard offline vertex cut \vtx, \ < 30 cm was required to select collisions well centered 
in the fiducial area of the central detectors and to avoid tracks close to the magnet poles. 

2. The maximum event multiplicity allowed was 15 tracks to suppress the contamination of 
non-UPC (mainly beam-gas and peripheral nuclear) reactions that fired the UPC trigger 

At variance with standard J/\\t ~> e+e^ analyses in nuclear AuH-Au reactions |40| which have to 
deal with large particle multiplicities, we did not need to apply very strict PID cuts in order to 
identify electrons in the clean UPC environment: 

1. RICH multiplicity «o >2 selects which fire 2 or more tubes separated by the nominal 
ring radius. 

2. Good CNT-EMCal matching is required for candidate tracks with an associated EMCal 
cluster without dead-warn towers within a 2 x 2 tile. 



3. An EMCal cluster energy cut (Ei > 1 GeV 1 1 £2 > 1 GeV) is applied to select candidate e 
in the plateau region above the turn-on curve of the ERT trigger (with 0.8 GeV threshold). 
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Beyond those global or single-track cuts, an additional "coherent" identification cut was applied 
by selecting only those e+e^ candidates detected in opposite arms (armi ^ arm2). Such a cut 
enhances the sample of back-to-back di-electrons with low pair pr as expected for yy, yA pro- 
duction. Finally, J/\\t were reconstructed by standard invariant mass analysis of all that passed 
the aforementioned analysis cuts. Any remaining background was removed from the m,„,, dis- 
tribution by directly subtracting the wrong sign pairs (e+e+ or e^e^) from the unlike-sign pairs 
in a bin-to-bin basis. The J/\\t yield is extracted directly by adding the number of counts within 
±3o (roughly 2.7 - 3.5 GeV/c^) of a Gaussian fit of the experimental mass peak. 



2.4 Acceptance and efficiency corrections 

Physical cross-sections are obtained after correcting the raw number of signal counts for: (i) 
the geometrical acceptance of our detector system, and (ii) the efficiency losses introduced by 
the aforementioned analysis cuts. Acceptance and efficiency corrections were obtained with a 
full Monte Carlo of the experimental apparatus with realistic input distributions of the physi- 
cal signals. We generated 10^ events for each one of the two coherent process of interest: J/\\f 
and high-mass e+e^ pairs production in AuH-Au collisions accompanied by Xn forward emis- 
sion, with the Starlight Monte Carlo 126112711411 . Such a model reproduces well the existing 
(f'N / dyddsfdpr distribution of coherent p production in UPC AuH-Au events measured at RHIC 
by STAR 1 28 1 . The coherent events were simulated in the PHENIX detector using GEANT 1 42 1 
and passed through the same reconstruction program as the real data. Figure|4]shows, for illustra- 
tive purposes, the invariant mass dN/dniee distributions of the coherent J/\\t's and e+e^ signals 
given by Starlight and a fast Monte Carlo version of our code |41 1. A Gaussian signal is ex- 
pected with experimental width of 1 00 Me V/c^ at m,>,v ~ 3 . 1 Ge V/c^ on top of a power-law-like 
e+e^ continuum. The correction factors obtained from our simulation studies are: 

1. Acceptance {J/\\t ^ e+e^ decay detected within \y\ < 0.5 and full (|)): Acc = 5.7%. 

2. Efficiency losses due to all cuts (including yield extraction within ±3a): e^™ = 56.4%. 

3. Unlike-sign background subtraction results in a negligible efficiency loss of the J/\\t signal. 



|T|(e'')| < 0.5, |r|(e")| <0.5 
Experimental J/'V width 




3.5 



4.5 



M,™ [GeV] 



Fig. 4. Expected invariant mass distribution of dielectron pairs from jA J/\\f e+e^ (red points) and 
YY^ e+e^ (histogram) as given by the Starlight model f41} for UPC Au+Au collisions at ^/s^ = 200 GeV 
detected in PHENIX. The hatched area indicates the region accessible with our ERT trigger. 



3 EXPERIMENTAL RESULTS 

Fig. |5] left, shows the invariant mass of unlike-sign pairs (red-filled histogram) and same-sign 
background pairs (yellow-filled histo) obtained from our analysis after application of the global. 
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single and pair cuts described in the previous section. The wrong-sign background is very small 
as expected from the MC results and has a flat pr distribution which extends far from the low-pr 
coherent region (Fig.|5]right). In what follows, all (m,„,, and dN/dpr) distributions have the like- 
sign background removed. Note that below m,„,, w 1.8 GeV/c^, the applied offline cuts above 
the ERT threshold energy (£1.2 > 1 GeV) are responsible for the effective lack of counts. 
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Fig. 5. Invariant mass (left) and pair pr (right) distributions of unlike-sign (red dots) and like-sign (open 
stars) e^e^ pairs identified in UPC Au+Au collisions at -v/i^ = 200 GeV. 



3.1 Invariant mass distribution 

The final measured background-subtracted t/A^/t/m^-g distribution is shown in Fig.|6lincluding the 
simulated e+e^ continuum power-law curve absolutely normalized at the measured dN/dniee 
yield at ntee = 2 GeV/c^ combined with a fit to a Gaussian at the J/\\f peak. The total number of 
counts for both physical signals is A^e±+i/y = 40 ± 6 (stat). Despite the poor statistics, the shape 
of the di-electron continuum is consistent with the power-law-like distribution expected for the 
reconstructed shape of the e+e^ MC signal (Fig. |4j. The dotted curves indicate the maximum 
and minimum e+e^ continuum contributions that we consider in this preliminary analysis and 
have been obtained by absolutely normalizing the power-law curve at the upper and lower values 
given by the ±lo yield uncertainty at the lUee = 2.0 GeV/c^ bin. From this plot it is obvious that 
the largest contribution to the systematic uncertainty in the extraction of the J/\\t signal comes 
from the subtraction of this physical background. 

Fig. shows the resulting invariant mass distribution obtained by subtracting the fitted power- 
law curve of the dielectron continuum from the total experimental e+e^ pairs distribution. There 
is a clear peak whose position and width are perfectly consistent with the expected signal 
from our full MC: nij/^ = 3.096 ± 130 MeV/c^ (the mass is right at the PDG value, though 
the width seems to be ^30 MeV/c^ larger than the simulated value). The total number of 7/i|/'s 
within ±3o of the peak position (2.7 - 3.5 GeV/c^) is: Nj/^ = 10± 3 (stat) ± 3 (syst.), where 
the systematic uncertainty is completely dominated by the di-electron continuum subtraction 
method. More detailed analyses are being currently carried out to get a better handle on the 
possible continuum contamination to the total J/\\t signal. 
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Fig. 6. Invariant mass distribution of e+e^ pairs measured in UPC Au+Au collisions at y^s^ = 200 GeV 
fitted to the combination of a dielectron continuum (power-law distribution) and a (Gaussian) signal. 
The two additional dashed curves indicate the maximum and minimum continuum contributions considered 
in this analysis below the y/\|/ signal region. 



3 

o 

e" - 

■o 

I 6- 



4- 



PHENIX AuAu UPC s^^ = 200 GeV 

• e*e" (coherent continuum subtracted) 



PHENIX Preliminary 



iPj,^ = 3.096 + 0.129 GeVIc' 

x'/NDF = 1.71/2.00 
N„ = 10±3(slat)±3(syst) 




5 6, 
iTiee (GeV/c ) 



Fig. 7. y/\|/ invariant mass distribution obtained subtracting from the total e+e^ pairs signal the fitted 
dielectron continuum shown in Fig.|6| The quoted number of J/\\l is just the number of counts within m^g 
= 2.7 - 3.5 GeV/c^ (the error is just statistical). 
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3.2 Pair transverse momentum distribution 

Fig.|8]shows the transverse momentum distribution of all the reconstructed e+e^ continuum pairs 
and J/\\t in UPC Au+Au collisions. Their spectrum is clearly peaked at low pj as expected for 
coherent production. The shape of the pr spectrum itself is well reproduced by a nuclear form 
factor fit l43l 



dNe, 

dp J 



' C ■ pt ■ \F {pt)\^ with F{pt) 



sin{R ■ pt) — R ■ pt ■ cos{R ■ pr) 
{R-pT?-{l+alp\) 



(10) 



with Au nuclear radius and diffusivity fixed to their known values, R = 6.38 fm aq = 0.54 fm 1441 . 
and with absolute normalization factor C as the only free parameter 




dN„/dp,= A*p,*|F(p,)|* 

sin(R'p^) - R'p/cos(RX) 
(R*p,f *(1+ajp?) 
Au distrib.: R = 6.38 fm, a. = 0.54 fm 



F(Pt) -- 



A = 1 300.0 ± 54.3 (z /NDF= 51 .9/32.0) 
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Fig. 8. Transverse momentum distribution of reconstructed e+e^ continuum pairs and y/\(/ from ultra- 
peripheral Au+Au collisions compared to the theoretical expectations of coherent photoproduction with a 
realistic nuclear form factor, Eq. <10t . 



3.3 Cross-section for coherent UPC production at mid-rapidity 

The final cross-section for coherent photoproduction at midrapidity in UPC Au+Au colli- 
sions at JT^ = 200 GeV accompanied by Au breakup is: 



d(3 



UPCJ/^ 



dy 



1 ^y/w 1 
= — = 48. ± 14. (Stat) ± 16. (syst) ^ib. 

v|<0.5 J^CC ■ Ereco ' ^trigg ' -^int '^3' 



where all correction factors (and corresponding uncertainties) have been obtained as described 
in previous sections, and BR = 5.93% is the known J/\\t di-electron branching ratio. The final 
cross-section is in very good agreement, within the (still large) experimental errors, with 
the theoretical values computed in |4r''221 as shown in Fig. |9] (the predictions of |22| have been 
scaled down by the nuclear breakup probability Pxn ^ 0.64). The current uncertainties preclude 
yet any detailed conclusion regarding the two crucial ingredients of the models (nuclear gluon 
distribution and J/\\t nuclear absorption cross-section). Whereas the systematical uncertainty of 
the measurement is linked to the (yet) imprecise contribution from coherent yy e+e^ pairs 
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below the peak, the statistical error can only be improved with a (much) higher luminosity 
Au+Au run. The contribution of the incoherent yN production to the total dOupc j/\^/dy\ ^ 

(included in the predictions of Strikman et al. 1221 but absent from Starlight L41J ') needs to be 
elucidated too. 
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9. Preliminary cross-section of coherent J/\\l production at mid-rapidity in UPC Au-l-Au collisions at 
= 200 GeV compared to the theoretical calculations 14 11221 . The error bar (box) shows the statistical 



(systematical) uncertainties of the measurement. 



4 CONCLUSION 

We have presented preliminary PHENIX results from an analysis of AuH-Au iy/s^ = 200 GeV) 
UPC triggered events aiming at the measurement of coherent photoproduction of J/\\t and high- 
mass e+e^ pairs in jA and yy processes respectively. Ultraperipheral AuH-Au collisions are 
tagged by forward (ZDC) neutron detection from the (single or double) Au* dissociation. Clear 
indications of J/\\i and high mass dielectron continuum have been found in the data. The total 
number of "physics" e+e^ pairs is: dN^±j_j j^/ dy = 40 ± 6 (stat). Their pj spectrum is peaked at 
low pt k, 90 MeV/c as expected for coherent photoproduction with a realistic nuclear form fac- 
tor After subtraction of the physical e+e^ signal, the measured invariant mass distribution has a 
clear peak at the mass with experimental width in good agreement with a full GEANT-based 
simulation for UPC production and reconstruction in PHENIX. The measured number of 7/i|/'s 
in PHENIX acceptance is: dN /dy = 10 ± 3 (stat) ± 3 (syst). After correcting for acceptance and 
efficiency losses and normalizing by the measured luminosity, we obtain a preliminary cross- 
section of d(3/dy\y=o = 48 ± 16 (stat) ±18 (syst) which is consistent, within errors, with 
theoretical expectations. Foreseen improvements in our analysis will likely reduce the experi- 
mental uncertainties of the measured cross-section and provide more quantitative information on 
the nuclear gluon distribution and absorption in cold nuclear matter at RHIC energies. 
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